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^ ■ Abstract 

H ■ 

We propose a method to compute the Isgur-Wise form factors Ti/2(1) and t^/2{^) 

for the decay of B mesons into orbitahy excited (P wave) D** charmed mesons on 
the lattice in the static hmit. We also present the result of an exploratory numer- 
ical simulation which shows that the signal/noise ratio allows for a more dedicated 
computation. We find Ti/2(1) = 0.38(5) and t^/2{^) = 0.53(8), with yet unknown 
systematic errors. These preliminary numbers agree fairly well with theoretical ex- 
pectation. 

PACS: 12.38.Gc (Lattice QCD calculations), 12.39. Hg (Heavy quark effective theory), 13.20. He (Lep- 
tonic/semileptonic decays of bottom mesons). 



Introduction 

The scalar heavy-light mesons and more generally the first orbital excitations D** have 
attracted attention since years and they still remain somehow mysterious. The recent 
discovery of a cs-scalar meson significantly lighter than expected has renewed the interest 
in these states [1,2]. There have been several lattice studies of this spectrum [3,4] and a 
recent rather complete one compares quenched and unquenched [5] computations. Recently 
the Hq -^ H-K transition (scalar-pseudoscalar-pion) have also been considered [6,7]. 

The transitions of the type B -^ D**lv raise a serious problem. In the infinite mass 
limit these decays are described by the Isgur-Wise form factors ri/2 and r3/2 [8]. To make 
a long story short, a series of sum rules [9]- [14] have been derived from QCD, all indicating 
that r3/2 should be significantly larger than ri/2- These sum rules relate the Tj form 
factors, as well as form factors related to excitations, to derivatives of the ground state 
Isgur-Wise function ^ and allow to bound the latter derivatives in an efficient and useful 
way [15]- [17]. Not only does the slope of ^ verify p^ > 3/4 but also the curvature and even 
higher derivatives are bound. The limit in which ri/2 = has been baptised "BPS" by 
Uraltsev [18]- [19] and was proven to provide intersting hints. 

However the theoretical prediction that r^ i > ^1^/2 ^-nd hence that the decay B ^^ D^ 
should be significantly larger than the B ^ Dq is not verified by experiment [2,20]. This 
is the '1/2 > 3/2' paradox [21]. One might incriminate the corrections to the infinite mass 
limit. Another possibility could be that the sum rules are fulfilled by higher excitations 
and that the ground state obeys an opposite hierarchy i.e. r^i < ''"1/2 • 

To answer to this question one needs to compute directly Tg i and rj; i- Here we propose 
a lattice method to do that. We will work in the static quark limit, tti^c -^ 00, with the 
four vectors v' = v = (1, 0, 0, 0), and we will exhibit we will exhibit operators whose matrix 
elements allow to measure these form factors. 

This letter is meant to propose this new method and to make a feasibility study. We 
do not intend at this stage to provide accurate results for these form factors but merely to 
describe the principle of the method and to show with preliminary simulations that there 
is good hope to make the precision calculation. 

1 Principle of the calculation 

We are concerned with the matrix element of an electroweak current between a pseu- 
doscalar or vector heavy-light meson H^*^ and an orbitally excited one H**. However, in 
the conditions of the infinite mass limit on the lattice with the heavy quarks at rest, both 
in the initial and final state (f^ = v'), this matrix element vanishes. 

The way out is to use a series of relations derived in ref. [22]. In that paper it has been 
shown that in the case of a matrix element which vanishes linearly in the difference v' — v, 
when v' -^ v, there are non-vanishing forward matrix elements (for v' = v) involving the 

^ There is no mathematical im.possibihty for the sum rules to be fulfilled with an reversed hierarchy for 
the ground state, but it does not seem very likely and is not seen in models. 



covariant derivative operator D^. These matrix elements are proportional to the infinite 
mass limit form factors ri (1) or rafl). 

2 2 

Let us summarise their proof using different notations: For simplicity we take v' = 
(1, 0, 0, 0), and v = v' + f_|_, where v_[_ is spatial up to higher orders in the difference v' — v. 
We assume that for some Dirac matrix Ti 

{H**{v')\h{v')Vih{v)\H^*\v)) = tTv^mT,{w) + • ■ • , (1) 

where w = v ■ v', j = 1/2, 3/2, and /,m = 1,3 are spatial indices, t'™ is a tensor which 
depends on the final state [H**) and the initial state {H* or H). The dots represent higher 
orders in v' — v. From translational invariance in the time direction, 

-ido{H**{v')\h{v')TMv)\H^*\v)) = (2) 



-i{H**{v')\h{v') 



T, D^ + D^ T, 



h{v)\H^*\v)) = tTv^^Tj{w){MH*,-MH) + 



The authors of ref. [22] use the field equation: {v ■ D)h{v) = 0, which implies that 

D^h{v') = 0, D%{v) = -{D-v^)h{v), (3) 

whence from (2) 

z{H**{v')\h{v')Ti{D ■ v^)h{v)\H^*\v)) = tTv^mT.iw) {Mh» - Mh) + ■ ■ ■ , (4) 

which has a finite limit when f _|_ ^ namely 

i{H**iv)\hiv)TiD"'hiv)\H^*\v)) = Cr,(l) (M^^ - Mh) . (5) 

Applying eq. (1) to the J = Hq state we get from ref. [8]: 

{H*iv')\A,\Hiv))^-T.{w)v^,, (6) 

where Ai is the axial current in the spatial direction i and where our states normalisation 
is l/v2M times the one used in ref. [8]. From eq. (6) it results that 

{H*,{v)\A,D,\Hiv)) =tg,, (M^. - M^) ri (1). (7) 

Analogously for the J = 2 H2 state we have 

{H;iv')\A,\H{v)) = V3Tiiw)e*'v^, + --- , (8) 

where e*- is the polarisation tensor, whence 

{H;{v)\AD,\H{v)) = -tV3{MH^ - Mh) Ml)e*^ . (9) 



2 Lattice calculations 

To compute the matrix elements in eqs.(7) and (9) on the lattice we first need a discretized 
expression for the covariant derivative. We choose the symmetrised form 



D,{x, t)^ — (Uiix, t) -Ul{x- 1 1) 



(10) 



where Ui{x,t) is the link variable of the lattice. 



2.1 Interpolating fields 

The interpolating fields for orbitally excited states have been studied in ref. [24] . Smearing 
is used not only to improve the signal/noise ratio by better isolating the ground state, but 
also to produce convenient interpolating fields for the 0~, 0"^ and 2^ states. Inspired by 
ref. [25] we replace the quark fields q{x) by 



q{x) 



+ 




x+(k-l)i) 



q{ 



x+ri) 



{ri) Sii 
1 



(11) 



where the upper (lower) expressions generate negative (positive) parity smearing functions. 
The vector (±r?) Sn is introduced to generate an orbital excitation in the direction /. The 
wave function 0(r) is a radial function chosen to optimise the overlap with the ground 
state. We take 0(r) = e~^'^'', where Rb is a parameter which is fixed by requiring the 
smearing to be optimal. Note that it is not necessary to normalize the wave function since 
the normalisation factors cancel in the computation of matrix elements. The smearing 
also includes the so-called fuzzing, see ref. [23]. For convenience, we will use the following 
notation for the interpolating fields: 



{h{x] 



lkr\x) 



q{x + f)} , 



(12) 



where r' indicates the presence of (±ri) Su in eq. (11), "t is a generic vector for the distance 
between the light and heavy quark field and F = 1 (F = 75) for the 0^ (0~) meson. 

Using the smeared quark fields from eq. (11) we now define the interpolating fields. We 
concentrate on the 0+ (2+) states which correspond to j = 1/2 (j = 3/2). The 0~''-state 
can be described according to two distinct interpolating fields: 



a) h{x) q{x + r) and 



b) —= h{x) (7 ■ r) q{x + r). 
V3 



(13) 



These two interpolating fields differ in that the Dirac matrix is diagonal (antidiagonal) for 
1 (7 ■ r) inducing the coupling of the heavy quark to the "small" ( "large" ) component of 



the light quark field. The latter is just the quark model combination of quark-spin 1 with 
orbital momentum 1 to generate J = O"*". 

Concerning the 2^ states, the same duality of interpolating fields exists. In this letter 
we only consider the quark-model type interpolating fields. This gives the five J = 2 
states [24] which we may write as follows: 



a) 



1 ^ 



\/2 



h{x) [7i ■ rj{x) + 7j ■ ri{x)] q{x + r), i ^ j, 



(14) 



b) -^h{x) [71 ■ ri{x) - 72 ■ r2{x)] q{x + r), 

v2 

1 _ 

c) -J^K^) [71 ■ n(a;) + 72 ■ r2{x) - 273 ■ r'i{x) ] q{x + r). 

These are, as expected, symmetric traceless tensors. 

2.2 Two-point Green functions and 1/2 — 3/2 mass splitting 

The 0^ two-point Green function is written as 



Cl 



2:0 



ETr 



,d 1 + 7o 



pu 



S{r{Q),Q-x + r{x)-U) 



(15) 



u 



when we use the interpolating field in eq. (13, a) ^. Pq^ is a temporal Wilson line ^ 
corresponding to the Eichten-Hill action for the static quark [26]: 



Pl,^=5{x~y) \[ul^^{x + tj), 



(16) 



^2 f^X 



using hypercubic blocking [27] - [29]. 

The two-point Green functions with 7jrj interpolating fields allow an interesting com- 
parison between the j = 1/2 and the j = 3/2 cases. They will contain terms of the general 
form 



Cfj\0,Q = (j^Tr 



lrrA^)PL 



1 + 70 



7fcr/(x)S'(f(0), Q]X + r{x); U) 



;i7) 



u 



where (i,j) = {k,l) for the case of interpolating field (14, a) ot i = j,k = I for the 
case (14,b), (14, c) and for the 0"^ case. J stands for the total angular momentum (J = 0, 2). 
After some simple algebra we can write 



-C^iO,Q = {J2T^PL 



{6,1 ± 2e^""a„) r,(0)^^rKx)5(f(0), 0; a; + r»; U) 



.at 



u 



^The superindex 1 (7 • r) refers to the use of the a (b) interpolating field in eq. (13). 
•^Indeed we compute the two point function using the interpolating fields in cqs. (13) and (14) properly 
shifted in space so as to have the light propagator ending at the origin. 



Let us define C2,5(r (o),r(x)) and C2,e'"(r (o),n(a;)) respectively the two terms in eq. (18). It 
is easy to see from tlie second interpolating field in eq. (13) that the O"*" two point Green 
function writes as 



'C'2;0 - 3 Z^ '^2,5(rj(0),rj(x)) + 3 Z^ [C'2,e'=(r,(0),r,(a;)) + C'2,efe(r,(0),r,(x))] 
i=l:3 i,j,k 

= C2,5{ri{0),ri{x)) + i [C'2,e3(ri(0),r2(a;)) " C'2,e3(r-2(0),ri(x))] , (19) 



where i,j, A; are in cyclic order and where we have taken advantage of the hypercubic 
symmetry in the r.h.s. 

Taking now any of the 2"*" meson interpolating fields and using the again the cubic 
symmetry we get 

■7 

—C2;2 = C'2,<S(ri(0),ri(x)) " " |_C'2,e3(ri(0),r2(a;)) — C'2,e»(r2(0),ri(x))J • (20) 

The difference between the j=l/2 and the j=3/2 state is thus related to the relative sign 
and coefficient of the e-term compared to the direct one. The effective energy is obtained 
by taking minus the time derivative of the logarithm of the two-point function. The energy 
difference between j = 1/2 and j = 3/2 is thus proportional to 



-I ( [ri(0)r2(x) - r2(0)ri(x)] TrPo° 



^3^^5(r(0), 0; X + f(x); U) 



(21) 
u 



In a non-relativistic limit f = ip/m. This imaginary velocity comes from the derivation 
versus the imaginary time. Then the expression in eq. (21) is reminiscent of a LS'-term: 
(rxp) -a except that the operators r and p are not taken at the same time, a in eq. (21) acts 
on the heavy quark but the trace will make it also act on the light quark. It is interesting 
that the coefficients of the last terms in eqs. (19) and (20) are in the ratio (1), (—1/2) which 
is exactly the ratio of the //^'-eigenvalues for j = 1/2, 3/2, built up from the combination 
of L = 1 and s = 1/2: 

2<L.S>=jij + l)-^-2 = ~^ for j = J/^ . (22) 

From eqs. (19) and (20) it is obvious that if these LS'-type terms did vanish the two-point 
correlators C2-o and C2;2 would be equal, which would then imply Mh* = Mh*. In this 
limit the normalisation of the interpolating fields in eqs. (13) and (14) has further ensured 
the equality of the multiplicative constants 22;0 and 22;2) where the Z2-j are defined from 

C2-AQ = i22;jfe-''"j'% (23) 

at large time tx- 



2.3 Three-point Green functions and ri-ra splitting 

The three point Green functions of the axial current using interpolating fields with '-fiVj 
will contain terms of the general form 

CtjA,Mty,t.) = UYlT^rl^^^M P^ ^7.75 (24) 

\ x,y 

UiiO, ty)Pl,^Sim, 0;x + f{x) + I; U) - U}i-l, ty)PiXSim, 0;x + f(x) - /; U) } 75' 



where we have used eq. (10) in units of a. Writing for short the term in the curly bracket 
as Di{y) ■ ■ -75, and we see can write 



-C^AJO, ty, t.) = Uj2T^T \{5,i ± ^e^^-a^) 75r,(0) P^^ ^^ A(y) ■ ■ ■ 75 



x,y 



(25) 



u 



where either {i,j) = {k,l) for the 0~ -^ 2+ transition, eq. (14, a), or i = j,k = I for 
the 0~ -^ 0~^ one and 0~ -^ 2+ with eqs. (14, b) and (14, c). Let us define C3^s{r{o),D{x)) 
and C3^em(r.(o),D;(a;)), respectively the two terms in eq. (25). From the interpolating field in 
eq. (13, b), using the fact that i = j and k = I, and choosing for simplicity k = 3, one can 
derive, 

-- 1 i 

-'^1,05 = -7^C'3,<5(r3(0),D3(x)) + -7= [C3,ei(r2iO),D3ix)) + C'3,€2(ri(0),D3{a;))] • (26) 

The axial matrix element is then given in Euclidean metric by 

< H*\A,D,\H >= ^2;o^2;5g3;05(0,t„t.) ^ ^^^^^ _ ^^^ ^^^^^ ^ ^27) 

'-'2;o(0, tj/) C2;5(0, tj; — tj^) " 2 

where we have used eq. (7). This leads, using cubic symmetry, to 

(M„, - Mu) ^/3 r. (1) = 2,„ Z,^ C,,(0,t,) C,,(0,t^ - t,) ' <'**' 

plus all terms deduced by cubic symmetry. 
From eq. (9), in Euclidean metric, we get 

" 2 ' C2;2(0,gC2;5(0,t.-y V2 ^ "^ ''' 2^^' ^ ^ 

where we have used the polarisation tensor 

71 
e=^ 0, (30) 





u 



for the 2^ state in eq. (14, a) with (?, j) = (1,2). From the interpolating field in eq. (14, a) 
using the fact that (i,j) = {k,l) = (1,2) one can derive, 
1 ^ 

-C3;25 = ^YTo ^ ^3,5{ri{0),Di{x)) + — /= [C3,e3(-n(0),D2ix)) + C'3,ea(r2(0),Di(x))] , (31) 

and using cubic symmetry and eq. (29) 

^ ' ' ^^ ^ ' ' C2;2(o,gc2;5(o,t,-g ^ ^ 

It can be checked that all the other states in eq. (17) lead to the same formula (32) up 
to a cubic rotation. The numerators in equations (28) and (32) exhibit identical C^^s 
terms and C^^e, the latter differing only by multiplicative coefficients which, once more, are 
proportional to the LS eigenvalues given in eq. (22). Combining the results of eqs. (19), 
(20), (28) and (32) we may conclude that, if the C2,e and C^^e terms did vanish, we would 
get Mh* = Mh* and T3/2 = ri/2. 

3 Condition and results of the simulation 

Results presented for Isgur-Wise functions ri (1) and rs (1) are obtained from the quenched 
simulation on a 16^ x 40 lattice at /3 = 6.0. We collected 580 independent SU{3) gauge 
configurations in the quenched approximation using the non perturbatively 0{a) improved 
Wilson fermion action with Csw = 1.769. The light-quark propagator is computed with 
the hopping parameter k, = 0.1334, which corresponds to a pseudoscalar "light-meson" 
mass of 800 MeV. For the static quark we use the "hyp" links as written previously. In 
fig. 1 we plot the binding energy for the scalar and the pseudoscalar meson and in fig. 2 
we plot the binding energy for the 2^ heavy-light meson. 

The scalar meson has been computed using the interpolating field ^ in eq. (13, a): 
h{x) q{x + r). The tensor meson has been computed using the properly averaged in- 
terpolating fields in eq. (14). We get A = tuh* - tuh = 400(12) - 411(16) MeV at 
(3 = 6.0. Only statistical errors are considered. It agrees reasonably with ref. [5] where 
A ~ 400(40) MeV. Our present signal for the tensor- meson effective mass is still very poor, 
rriH* — iTiH = 0.50(8) GeV, which leads to rriH* — thh* = 0.10(8) GeV. The large relative 
error reflects the poor quality of the plateau in flg. 2. Clearly a more reflned simulation is 
needed here. In particular we have not yet optimised the wave function for the smearing 
of the tensor meson. Our result agrees with the result of [5] where we read from table 2 
(Qs) : rnn* - thh = 0.48(2) GeV, and uih* - mn* = 0.08(4) GeV. 

Experimentally the situation is not yet clear: whereas Belle [2] reported mo* — mD* = 
153(36) MeV ^ FOCUS [30] flnds tud* - mo* = 61(41) MeV. Anyway large ifrric correc- 
tions are expected. 

"^This choice shows up a better signal than the one using eq. (13, b). A comparison of these signals has 
been performed in [7]. 

^Note however than in Belle experiment the narrow and broad J^ = 1+ resonances, usually interpreted 
as j — 3/2, j — 1/2 respectively, are practically degenerate in mass: m^jo = 2421(2) MeV and nijj'o = 
2427(50) 
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Figure 1: Signals for the effective binding energies for the pseudoscalar and the scalar heavy-light 
mesons. 




Figure 2: Signal for the effective binding energies E^s for the 2^ heavy-light mesons. 



In fig 3, we plot the ratios 
n(l) = 

2 ^ ' 



1 



^2;0 ^2;5 C'3;05(0, tj^, t^ 



(Mh* - Mn) C2;o(0, tj,) C2;5(0, t. - ty) ' 

2 1 22;2 ^2;5 C*3;25(0, ty, ta;) 



(33) 



3 (Mh* - Mh) C2;2(0, tj,) C2;5(0, t. - ty) ' 

where the source operator has been fixed at t^ = 13a. The equality is valid on the plateau. 
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Figure 3: Signals for the ratios defined in eq. (33); from the fit in t/a G [4,9] and t/a G [3,5] 



respectively, we obtain the value of t\ and ts . 

2 2 



4 Results, Discussion and Conclusions 

In this letter we propose a method to compute on the lattice, at the infinite mass limit, the 
zero recoil Isgur-Wise form factors ti(1) and T3{1) relevant to the decay of a heavy pseu- 
doscalar meson into orbitally excited states. The main feature of the method is contained 
in eqs. (7) and (9). It uses matrix elements of the axial current multiplied by covariant 
derivatives. 

We have also performed an exploratory lattice study in order to estimate if this method 
is practically usable. We find that the signal/noise ratio is encouraging if one considers 
that there is still room for improvement. 
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Our results are 

ri(l) = 0.38(4)(?) and r|(l) = 0.53(8)(?) , (34) 

where the question mark represent yet unknown systematic errors. We have also: 

r|(lf-r|(l)^ = 0.13(8)(?) (35) 

Within 1.5(7 it saturates the Uraltsev sum rule [10]: Yin \'^i V-)? ~ ki (1)P = i- Note 

2 2 

that an approximate saturation of the sum rule by the ground states is seen in several 
models [31,33] although there is no strong theoretical reason for that. 

The result for ti(1) = 0.38(4)(?) is presumably more reliable than the one on r3(l) 
since the two-point signal for the O'^ meson is much better than the one for the 2+ meson, 
see figs. 1 and 2. Our result for ri(l) is somewhat larger than the predictions of the 
covariant quark models a la Bakamjian-Thomas (BT) [31] which predict, for the preferred 
potentials ^ ri(l) G [0.1,0.23] andr3(l) G [0.43,0.54]. The latter agrees well with eq. (34). 
Both numbers of eq. (34) are compatible with a recent calculation based on a covariant 
light-front approach with simple harmonic oscillator wave functions which are not derived 
from a potential ^: ti (1) = 0.31 and T3(1) = 0.61 [33]. It is interesting to not that both in 
ref. [31] and in [33], ts (1) — ri (1) ~ 0.3, which might be a general feature of BT covariant 
quark models (see eq. (5.1) in ref. [31]). This is somewhat larger than the difference 
between central values of (34). We also agree with an older QCD sum rule estimate [32]: 
ri(l) = 0.35(8). 

The quantity we compute on the lattice is a physical quantity: the product of a mass 
difference times the form factors ti(1) and T3(1). We thus expect no multiplicative renor- 
malisation to be needed. A closer scrutiny of this question is underway, in particular to 
understand if the hypercubic treatment of the Wilson line can have some effect on the 
discretized covariant derivative we use. Of course a complete control of systematic effects 
is also needed: finite volume, mass of the light quark, finite lattice spacing. In section 2.1 a 
duality of interpolating fields has been pointed out. A systematic comparison of the their 
predictions is still missing. 
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®It should be stressed that the BT method provides a framework in which different potentials can be 
used. The physics prediction depends, of course, on the chosen potential. 

^The covariant light-front framework is equivalent, in the infinite mass limit, to the BT one. The 
practical predictions depend, however, on the chosen parameters and shape of the wave function. We 
worry if the use of gaussian wave functions, which is frequent, is a good one as it neglects the short 
distance potential. 
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